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Abstract. To set useful limits on the abundance of small-scale dark matter 
halos (subhalos) in a galaxy scale, we have carried out mid-infrared imaging 
and integral-field spectroscopy for a sample of quadruple lens systems showing 
anomalous flux ratios. These observations using Subaru have been successful 
for distinguishing millilensing by subhalos from microlensing by stars. Current 
status for our lensing analysis of dark matter substructure is reported. 



1. Introduction 

Recent high-resolution N-body simulations based on Cold Dark Matter (CDM) 
theory highlight a so-called missing satellite problem, i.e., CDM predicts the 
existence of more than several hundred dark satellites (or subhalos) in a galaxy- 
sized halo, in sharp contrast to the observed number of about 20 Milky Way 
satellites (e.g., Diemand et al. 2007 for recent studies). To clarify this issue, 
gravitational lensing offers us an invaluable insight into such numerous CDM 
subhalos that reside in lensing galaxies. In particular, anomalous flux ratios in 
lensed QSOs, namely those hardly reproduced by any lens models with a smooth 
density distribution, are of special interest, since lens substructures are able to 
cause such flux anomalies (e.g., Metcalf & Madau 2001; Chiba 2002). 

Here, we report on our Subaru observations of a sample of quadruple lenses 
with anomalous flux ratios, based on mid-infrared imaging and integral- field 
spectroscopy. Mid-infrared imaging of lenses is advantageous because the flux 
is free from differential extinction among different images and it is also free 
from microlensing by stars. Integral-field spectroscopy provides both spatial and 
spectral information simultaneously on each lensed image. We select eight QSOs 
with four lensed images, PG1115-h080, B1422-F231, MG0414-h0534, Q2237-h030, 
H1413-hll7, HS0810-h2554, and WFI2026-4536 for mid-infrared imaging, and 
RXJ1131— 1231 for integral-field spectroscopy. We briefly report the results for 
PG1115-F080 and 61422-^231 (Chiba et al. 2005) and RXJ1131-1231 (Sugai 
et al. 2007); those for other targets will be published elsewhere (Minezaki et al. 
2008 in preparation). We present new limits on substructure in the observed 
lensed systems and implications for a missing satellite problem are discussed. 
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2. Targets and Observations 

PG1115+080 at redshift zs = 1.72 and B1422+231 at zs = 3.62 are lensed 
by foreground ellipticals at = 0.31 and 0.34, respectively. The former lens 
system holds the closely separated pair of images Al and A2 with a separation 
of 0."48, and this configuration emerges if the QSO is close to and inside a fold 
caustic provided by the lens (Figure 1). The latter shows the colinear, three 
highly magnified images, A, B, and C, and this configuration emerges if the 
QSO is close to and inside a cusp caustic. In such lens systems associated with 
a fold or cusp caustic, there exists a universal relation between the image fluxes, 
i.e., A2/A1= 1 or (A+C)/B= 1, whereas the observed optical flux ratios violate 
these rules significantly, A2/A1= 0.64 ± 0.02 and (A+C)/B= 1.50 ± 0.01. 

The mid-infrared imaging of PG1115-I-080 and B1422-I-231 was carried out 
on the nights of UT 2004 May 5 and 6, using the cooled mid-infrared camera 
and spectrometer (COMICS). The field of view is 42" x 32" and the pixel scale is 
0."129 pixel"^. We used the Nil. 7 filter, whose effective wavelength and band- 
width are Ac = 11.67 fim and AA = 1.05 iJ.m, respectively. The FWHM of PSF 
was 0."33 at small airmass. The observed mid-infrared waveband corresponds 
to the near-infrared waveband in the rest frame, and its flux is dominated by 
thermal radiation from hot dust located at the innermost region of a dust torus. 
The inner radius of a dust torus, which is determined by the highest sublima- 
tion temperature of dust (T ~ 1800 K) and the UV luminosity of a QSO central 
engine, is generally much larger than Einstein radii of foreground stars, so that 
the observed mid-infrared flux is free from microlensing effects. 

RXJ1131— 1231 is unique in its low redshift of a source image zs = 0.658 
lensed by an elliptical at = 0.295. The lens shows three roughly co-linear 
images. A, B, and, C, being characteristic of a cusp singularity, where A is 
brightest (Figure 2). The observed flux ratios, (B-|-C)/A~ 2.1 in the V band 
and 2.2 in the R band, deviate significantly from the rule (B-|-C)/A= 1. 

Using the IFS mode of the Kyoto tridimensional spectrograph II (Kyoto 
3DII), we observed RXJ1131-1231 on the night of UT 2005 February 8. The 
IFS mode uses an array of 37 x 37 lenslets, enabling us to obtain spectra of ~ 10^ 
spatial elements. The spectral range from 7300A to 9150A was observed in each 
of two one-hour exposures. With the spatial sampling of 0".096 lenslet^^, the 
field of view of ~ 3" covered the three bright lensed images. We measured the 
emission-line fluxes of both the BLR H/3 and the NLR [OIII]AA4959,5007 for 
images A, B, and C simultaneously. The H/3 and [OIII] lines are very close in 
wavelength, so that the effect of differential reddening between them is negligible. 



3. Results 

The mid-infrared images at A = 11.7 fim of PG1115-f080 and B1422+231 are 
presented in Figure 1. It is evident that the lensed images in concern, (Al, A2) 
and (A, B, C), respectively, are clearly detected and well separated from each 
other. The mid- infrared flux ratios are given as A2/Al= 0.93 it 0.06, B/Al= 
0.16±0.07, and C/A1= 0.21±0.04 for PG1115+080, and (A+C)/B= 1.51±0.06, 
A/B= 0.94 ± 0.05, and C/B= 0.57 ± 0.06 for B1422-h231. This suggests that 
the mid-infrared flux ratio for A2/A1 of PGl 115-1-080 is consistent with the 
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Figure 1. Left: Mid-infrared images at 11.7 /xm for PGlllS+OSO (top) and 
B1422-I-231 (bottom). The north is up and the east is left. Right: Smooth 

lens models for these systems, where solid and dotted lines denote the caustics 
and critical curves, respectively, and filled circles indicate the source positions. 



prediction of a smooth lens model (~ 1), while the optical flux ratio is much 
smaller (~ 0.65), and that the mid-infrared flux ratios for (A, B, C) images of 
B1422-I-231 remain anomalous, contrary to the prediction of a smooth lens. 

Based on the dust reverberation method for estimating the size of a dust 
torus (Minezaki et al. 2004), we obtain an angular size of a source image as 
6's ^ 1 X 10"^ arcsec and 3.7 x 10"'' arcsec for PG1115+080 and B1422+231, 
respectively. For PG1115-I-080, any substructure causing its optical anomalous 
flux ratio A2/A1 should have a small Einstein angle 9e compared to 9s, because 
its mid-infrared ratio remains unaffected. This suggests that a substructure mass 
inside denoted as Me, should be smaller than 20 M(:j, being comparable to 
the mass of a star, i.e., microlensing causes an optical anomalous flux ratio. 
For B1422-I-231, the presence of flux anomaly even in mid-infrared waveband 
suggests Me ^ 200 Mq, i.e., lensing by a subhalo is most likely. 

Figure 2 shows the Kyoto 3DII spectra of RXJ1131 — 1231 for three bright 
images A, B, and C. It is interesting to remark that images B and C in [OIII] 
show nearly the same fluxes relative to image A, as expected from a smooth 
model at its cusp singularity. Thus, the absence of substructure lensing effects 
on this NLR [OIII] sets important limits on the mass of any substructures along 
the line of sight, as Mg < 10^ Mq. In contrast, the H/? line emission, which 
originates from the BLR, shows an anomaly in the flux ratio between images 
B and C, i.e., a factor two smaller C/B ratio than predicted by smooth-lens 
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Figure 2. Left: Kyoto 3DII spectra of RXJ1131-1231 for image A (solid 
line), B (dashed line), and C (dot-dashed line), where each has been ex- 
tracted with an 8-lenslet, i.e., a circular aperture with a diameter 0".77. Right: 
Smooth lens model for this lens system. 

models. The ratio of A/B in the H/3 line is well reproduced. The anomalous 
C/B ratio for the H/3 Une is caused most likely by microlensing of image C, with 
Me ^0.1 Mq for the mass of a substructure near image C. Wc have also found 
the slight difference of the H/3 line profile in image A from those in the other 
images, which suggests the presence of a small microlensing effect on image A. 

4. Prospects 

In addition to the above lens systems, we have already observed, using COMICS, 
MG0414+0534, Q2237+030, H1413+117, and HS0810+2554, and the calibra- 
tion and analysis arc underway. Our final target for mid-infrared imaging of 
quadruple lenses, WFI2026— 4536, is scheduled for observation using Gemini- 
South this year. Our preliminary statistical model implies that about 30 % of 
the lens systems (2 — 3 out of 8) show a flux anomaly with > 25 %, if subha- 
los are modeled by tidally-truncated singular isothermal spheres with a mass 
function predicted by the N-body simulations. Once the concrete observational 
information for all of these targets is ready, we will be able to set a more reliable 
constraint on the abundance of CDM subhalos. 
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